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We study the gravitational instability of an isothermal gaseous slab formed by cloud- 
cloud collision and compression at the cloud interface. The compressed gaseous slab 



p . rotates and has velocity shear except when the collision is not exactly head-on. The 

JL I effects of the rotation and velocity shear on the gravitational instability have been eval- 

5-1 I uated for the first time. We obtained the growth rate of perturbations as a function of 

^ ■ the wavelength for various gaseous slab models. We also obtained the changes in the 

density and velocity due to the perturbations. Two types of unstable modes are found; 
one is due to the self-gravity of the gaseous slab and the other is due to the velocity 
/\ • shear. The former instability leads to the formation of star and star clusters. Velocity 

j^ ■ shear decreases the growth rate of the former instability while it increases that of the 

latter instability. Rotation also decreases the growth of the former instability. As the 
slab grows in mass, the growth rate of the self-gravitational instability increases and 
eventually becomes greater than the growing timescale of the slab. Then a gravitation- 
ally bound condensation is formed and its mass is larger when the slab has stronger 
velocity shear and rotation. 
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1. Introduction 

The high- velocity cloud-cloud collision has been 
thought to be one of the important triggering 
mechanisms of the star and star cluster forma- 
tion. When two clouds collide with supersonic 
velocities, dense gaseous slab bounded by shock 
waves forms at the interface of the colliding 
clouds. The dense gaseous slab grows in mass 
and becomes unstable against fragmentation be- 
cause of the self-gravity. The fragments collapse 
and evolve into stars and star clusters. 

Many evidences have been found for star for- 
mation triggered by the cloud-cloud collision. 
NGC1333 (Loren 1976), GllO-13 (Odenwald et 
al. 1992), W49N (Serabyn et al. 1993), and 
Sagittarius B2 (Hasegawa et al. 1993) are nom- 
inated as the present sites of massive star for- 
mation triggered by cloud collision. The num- 
ber density of Galactic H II regions is propor- 
tional to the square of the local gas density, pu2 ■ 
This static suggests that substantial fraction of 
OB stars exciting the H II regions are formed by 
cloud-cloud collision (Scoville et al. 1986). 

It is also proposed that the populous and glob- 
ular star clusters are formed by high- velocity gas 
clouds collisions. Fujimoto and Noguchi (1990), 
Fujimoto and Kumai (1991a, b), Ashman and 
Zepf (1992), and Kumai, Basu and Fujimoto 
(1993) assessed the gas-dynamical circumstances 
of the Magellanic Clouds, M31, M33, M82, and 
other galaxies having young populous and glob- 
ular star clusters, and found that these galax- 
ies have large amount of gas in large-scale un- 
organized motion induced internally or exter- 
nally. In a great contrast to these galaxies, our 
Galaxy has neither young globular clusters, nor 
large-scale unorganized motion of gas. From 
this comparison, they developed the model that 
unorganized motion with velocity greater than 
~ 50km s~^ creates strongly-compressed regions 
of gas through cloud-cloud collision, and trig- 
ger the star cluster formation (see, Gunn 1980, 
Sabano and Tosa 1985, and Kang et al. 1990 



for another model of star cluster formation by 
cloud-cloud collision). In this model, the gas in 
the compressed slab fragments into clumps whose 
mass is comparable to that of the whole star clus- 
ter, and they cascade to smaller and smaller frag- 
ments, finally they form a gravitationally bound 
stellar systems, - globular clusters. 

In the star and star cluster formation model by 
cloud collision one must take account of oblique 
collisions since they are much more frequent 
than head-on collisions. When two clouds col- 
lide obliquely, the shock compressed gas slab ro- 
tates and is stressed by the velocity difference in 
the tangential component (velocity shear). The 
rotation and velocity shear may affect the frag- 
mentation of the shock compressed gas slab. The 
present paper evaluates the effects on the gravi- 
tational instability of the slab for the first time. 
It is shown that velocity shear increases the effec- 
tive sound speed of the gas and accordingly the 
minimum mass for instability (the Jeans mass). 
Rotation also increases the minimum mass for 
the instability. 

In section 2 we describe our slab-model which 
mimics the shock compressed gas formed by 
oblique collision of clouds. Our slab takes ac- 
count of rotation of the slab and velocity shear 
in the slab. In section 3, we formulate a linear 
analysis of perturbation around the equilibrium 
state. The numerical results are shown in sec- 
tion 4 consisting of subsections 4.1 through 4.4 
each of which are for non-rotating slab without 
shear, non-rotating slab with shear, rotating slab 
without shear, and rotating slab with shear, re- 
spectively. In section 5 we estimate the Jeans 
mass to apply the dispersion relation obtained in 
the previous section to star and star cluster for- 
mation. Section 6 is devoted to discussion and 
conclusion. 



2. Dynamical Models for a Shock-Com- 
pressed Slab of Self-Gravitating Gas 

In the supersonic cloud-cloud collision, a shock- 
compressed gaseous slab is formed at the inter- 
face of the two colliding clouds, (e.g., Stone 1970; 
Struck-Marcell 1982; Sabano and Tosa 1985). 
Figure 1 is the snapshot of a model of the oblique 
collision, with the Cartesian coordinates referred 
to the gaseous slab, where x — y plane is taken on 
the interface and the z-axis is taken normal to it. 
The gaseous slab rotates around the y-axis with 
angular velocity fi, and the compressed gas flows 
parallel to the x axis with velocities uq directed 
toward x = -|-oo in the space z > and x = — oo 
in z < 0. The shear flow is assumed to be uni- 
form in the x- and y- directions for simplicity, 
uq = uo{z). As shown in (c) of figure 1, the flow 
parallel to the x-axis is dominant over that to 
the z-axis, because the latter component is much 
reduced by the strong compression. We neglect 
the latter component in the present paper. 

In order to simplify the problem, we take the 
following three approximations. First we neglect 
the high-temperature region immediately behind 
the shock wave and assume the whole shock- 
compressed slab to be isothermal. The high tem- 
perature region is small in mass because of fast 
radiative cooling therefrom. Second the gaseous 
slab is assumed to be uniform in the x- and y- 
directions. As noticed by Stone (1970), the loss 
of gas across the lateral edges is negligibly small 
during the period of the collision. Third the 
gaseous slab is assumed to be in equilibrium, al- 
though its column mass actually increases with 
time. This equilibrium state can be achieved, 
when the sound wave propagates from the shock 
front to the interface before the gaseous slab 
grows considerably. As Welter (1982) showed, 
it is actually realized in the supersonic collision 
of gas clouds. When half of the slab thickness 
and isothermal sound speed are denoted with L 
and c respectively, we have for the sound crossing 



Fig. 1. — The obliquely colliding clouds, a: 
Snapshot of the ongoing collision. Thick arrows 
are a reference frame. The clouds move along 
horizontal lines and collide. We find a shock com- 
pressed slab at the interface of two clouds, b: A 
model of the shock-compressed slab (dark patch) . 
Circles spread by a light patch denote the pro- 
jected clouds on the x — z plane, c: Gas flows in 
and out of the slab. The pre-shocked gas streams 
along the horizontal lines, and compressed by the 
oblique shock. Although the post-shocked gas 
flows along dashed arrows, we approximate the 
stream line by solid arrows parallel to the slab. 



time t, between < z < L, 
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Similarly the increasing time tm of the slab mass 
is written as 



tr 
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M^t, , (2.2) 



where a, po, pi, and M_l denote respectively 
the column mass of the slab, preshock density, 
postshock density, and the Mach number of the 
gas incident on the shock front normally. Here 
we have used relations, a ~ 2piL, da/dt ~ 
2pQcM±, and pi/po = M±'^. We have thus quasi- 
equilibrium when tg < tnu and Mj_ > 1, both of 
which hold in the present supersonic collision of 
gas clouds. 

3. Normal Mode Analysis for a Plane- 
Parallel Compressed Gaseous Slab in 
Rotation and Nonuniform Parallel 
Streaming 

We analyze linear stability of a self-gravitating 
gaseous slab with the shear along the x-axis and 
the rigid rotation around the y-axis (figure 1). 
The basic equations governing the motion of gas 
are written as. 
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+ V • (pu) = 0, 


(3.2 


AV' = 4ttGp, 


(3.3 



where u,fl,r,ip,p,p are usual symbols, repre- 
senting the velocity of gas, angular velocity of 
the slab, position vector, gravitational potential, 
pressure, and gas density, respectively. As noted 
in section 2 and shown in figure 1, the angular 



velocity of the slab is parallel to the y-axis, and 
we have explicitly, 



n = {o,n,o). 



(3.4) 



Since we assume an isothermal gas, the pressure 
is written as 



p = c p. 



(3.5) 



The dynamical quantities u, p, tp,p are expressed 
as the sum of those in an equilibrium state and 
the deviations from it, the former and latter of 
which are specified with suffix and with hat and 
suffix 1, respectively, like for example, 

p{r,t)=po{z)+p^{r,t). (3.6) 

The velocity in equilibrium, Uq, has only the x- 
component and is a function of z, 



Wo = [uo{z), 0, 0] 



(3.7) 



Note that the other equilibrium quantities po and 
tpo are also functions of z and independent of t, 
X and y, since we have assumed that the slab is 
homogeneous in the x- and y-directions. 

3.1. Equilibrium Configuration of the 
Rotating Gaseous Slab vi^ith Shear 

Substituting equation (3.6) and similar equa- 
tions for other variables, into equations (3.1), and 
(3.3), and using equations (3.4), (3.5), and (3.7), 
we derive the equations for the equilibrium state, 
by omitting the terms with hatted symbols, 
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dz dz 



and 
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(3.8) 



(3.9) 



Equation of mass conservation (3.2) is automati- 
cally satisfied. Equation (3.8) represents that the 
Coriolis force per unit mass balances the gravita- 
tional acceleration and pressure gradient. Three 
variables, uo,po, and ipo are constrained by only 
two equations (3.8) and (3.9), and hence we have 



infinitely large number of solutions. In this paper 
we restrict ourselves to the case that the density 
distribution has the same functional form to that 
given by Spitzer (1942), 



po{z) = pcsech^ [jr) 



Then we obtain, 



uo{z) = [/-^tanh ij^j, 
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with nUH/c^ < 1, and 
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(3.12) 

(3.13) 
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where pc is the density at z = 0, and Hr is a 
density scale height which tend to i:f as i7 to 
zero. The velocity shear is duo/dz = U/H on 
the z = plane. The column density a is given 

by 

a = j pdz = 2pc^_Rtanh f — — | . (3.15) 

3.2. Linear Perturbed Equations 

Linearized equations of equations (3.1) to (3.3) 
and (3.5) are written as, by using equations (3.4), 
(3.6) and (3.7), 

duxi , dua:i . ^^0 , oo- 
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+Tr (pouzi) = 0, 
dz 


(3.19) 
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Since the equilibrium quantities are functions 
only of z, the perturbed ones can be expressed 
in the form of, for example, 

pi{r,t) = pi{z)e-x.Y)[i{kxX + kyy - ujt)] . (3.22) 

We introduce new variable y{z), consisting of 
four component, yi through 2/4 
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where gzi denotes the perturbed gravity in the 
z-direction, 

- - 1 ■ '-^) 

After straightforward manipulation we have for 
equations (3.16) through (3.21), 
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where the 4x4 matrix A has following elements 
1 dpo 2Qkx 
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(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 
^32 = 



and others are zero, i.e., A24 
^33 = ^42 = ^44 = 0. 

We integrate equation (3.25) and seek so- 
lutions satisfying the boundary conditions de- 
scribed in subsection 3.3. Then we obtain the 
dispersion relation, i.e., a; as a function of kx and 
ky for a given slab model. 

3.3. Boundary Conditions 

Since the gaseous slab is sandwiched and 
pressed together with two plane-parallel shock 
fronts, the boundary surface of our integration 
is assumed to be rigid and fixed [e.g. Stone 
(1970), Voit (1988)]. That is, we take into ac- 
count the fact that the plain shock front is stable 
against perturbation (e.g. Landau and Lifsitz 
1987). When a wavy perturbation is imposed 
along the shock front, the velocity component 
of gas parallel to the disturbed surface does not 
change across the shock front, but, the compo- 
nent normal to it decreases considerably due to 
the compression of gas. The resultant motion of 
the gas is divergent from the convex region (seen 
from the downstream) and similarly convergent 
to the concave one, restoring the distorted front 
surface to its original plane surface. 

We allow of the pressure discontinuity 6p due 
to the perturbation. According to equation (3.23), 



yi{L) =6p. 



(3.36) 



The perturbed velocity normal to the surface, Uzi 
or 1/2, must be zero. 



y2iL) = 0. 



(3.37) 



The perturbed potential ^1 is continuous and 
smooth at z = ibL. Solving the Poisson's equa- 
tion in the space outside the slab, A-f/'i = 0, we 
have 

Mz)=ML)exp[-^kx^ + ky\z-L)], (3.38) 

in z > L, and thus 



^4^,{L) = -^fkJTky'ML). (3.39) 

Using equations (3.23) and (3.24), we have the 
boundary conditions on 1/3 and 2/4, 



ysiL) = pc^i{L), 



(3.40) 



and 



yA{L) = -pc^kx^ + ky^ML). (3.41) 



Equations (3.36), (3.37), (3.40), and (3.41) form 
a complete set of the boundary conditions at 
z = L. Note, however, that it includes two pa- 
rameters, 5p and ^i(L), whose numerical choices 
are still free. 

Similar boundary conditions are imposed at 
z = —L. In total the boundary conditions have 
four degrees of freedom as well as the differential 
equation (3.25) has. 

3.4. Numerical Methods 

Our method to obtain the numerical disper- 
sion relation is essentially the same as those of 
Nakamura et al. (1991) and Matsumoto et al. 
(1994). We seek uj for which a solution of equa- 
tion (3.25) satisfys the boundary conditions at 
z = ibL using a bisection method. 



The general solution satisfying the boundary 
condition at 2; = L is expressed as a linear com- 
bination of two linearly independent solutions, 
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where y^'^^z) and y^'^'{z) are solutions of equa- 
tion (3.25) and their values at 2; = L are given 

by 

2/«(L) = 
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and 
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(3.44) 
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respectively. By the Runge-Kutta method we in- 
tegrated equation (3.25) to obtain y^'^^z) and 



y 



mi 



The parameters, 5p and ^{L) are left 
unfixed at this moment. 

Similarly we obtain the general solution, 
y{z) = 5p'y^^>{z) + ilj{—L)y^^>{z), which satis- 
fys the boundary condition at z = —L. The two 
general solutions coincide each other at z = 0, 
and thus 



cii 
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with 



yi^^\0) i,i = 1,2,3,4. 



(3.46) 



Equation (3.45) has a non-trivial solution when 
and only when the determinant of the matrix 
(cjj) vanishes. 



det 
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(3.47) 



We computed the determinant as a function of 
CO and sought lo for which the determinant has 
a sufficiently small value. See Matsumoto et al. 
(1994) for the algorithm to obtain cv for the de- 
terminant has a smaller value successively. 

4. Gravitational Instability of the Plane- 
Parallel Gaseous Slab 

The dispersion relations between uj and k[= 
{kx, ky)] are obtained for four cases with/without 
the shear and the rotation. First we study 
the gravitational instability of non-rotating slab 
without velocity shear in subsection 4.1. The ef- 
fects of velocity shear and rotation are shown sep- 
arately in subsections 4.2 and 4.3, respectively. 
We discuss the gravitational instability of a ro- 
tating slab with velocity shear in subsection 4.4. 

We represent all the dispersion relation in a 
nondimensional form for extensive applications 
of the results. For this purpose we normalize 
velocity, time, and density, with the isothermal 
sound speed c, the free-fall time l/y/27rGpc, and 
the density pc at z = 0, respectively. Thereby 
the length is automatically normalized with the 
density scalehight, H = c/^y2'irGpc- After the 
normalization, parameters are reduced to k, L, U, 
and fi. 

4.1. Non-rotating Slab without Velocity 
Shear (n = U = 0) 

Figure 2 shows the squared growth rates (—1^^) 
of the unstable perturbation as functions of kx 
for ky = 0. The thin and thick curves de- 
note the dispersion relations for L = 0.5 and 
1.0, respectively. The growth rate vanishes at 
kx = and the perturbation is stable when the 
wavenumber exceeds the critical one, kx > k^r- 
The critical wavenumber is kcr = 0.74 and 0.95 
for L = 0.5 and 1.0, respectively. The growth 



rate has its maximum —tOn 



0.17 and 0.34 at 



kxma.x = 0.37 and 0.46, respectively for L = 0.5 
and 1.0. 

Figure 2 suggests that — Wmax^ and fc^max de- 



crease as the slab thickness 2L decreases. Fig- 
ure 3 confirms it. For a given kx the mass con- 
tained in a wavelength decreases and hence the 
growth of the gravitational instability weakens as 
the slab becomes thinner. When the slab is thin, 
< L < 0.5, the dispersion relation shown in fig- 
ure 2 and 3 can be reproduced approximately by 
the gravitational instability of the infinitely thin 
sheet of gas whose column density is given by 
equation (3.15) (e.g. Stone 1970, and Goldreich 
and Tremaine 1979): 
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(4.1) 



Fig. 2. — Dispersion relations of the gravitational 
instability when U = i^ = 0. Thick and thin 
lines are results for the fixed and free boundaries, 
respectively. We write adopted L in the figure. 



Fig. 3. — The dependence of the maximum 
growth rates and their wavenumbers on the 
thickness of the slab for U = Q = ky = 0. The 
meanings of the thick and thin lines are the same 
as in figure 2. Upper two lines represents squared 
maximum growth rates, while lower ones are the 
wavenumbers at the maximum growth rates. 



with 



and 
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(4.3) 



[We suspend the normalized expression for quan- 
tities in equations (4.1) through (4.3).]. Also for 
a larger L equations (4.2) and (4.3) reproduce the 
qualitative features of the dispersion relation but 
less accurately. At the limit of L ^ oo, — Wmax^ 
and kxma.x tend to 0.45 and 0.47, respectively 
(See also Simon 1965), while equations (4.2) and 

(4.3) give -Wmax^ = 1-0 and /c^rmax = 1-0. 

The top diagram of figure 4 represents a cross 
sectional view of the fastest growing density per- 
turbation with the growth rate, — Wmax^ = 0.17, 
for the slab of L = 0.5. It delineates the one 
wave-length region -n/kxruax < x < Tr/kxma.x- 
In order to see the vertical structure, we magnify 
the slab thickness about ten times in the bot- 
tom diagram. Contours denote the gas density. 
Arrows in the bottom panel shows the velocity 
distribution. We find that the perturbed veloci- 
ties are predominantly parallel to the x-axis, and 
negligibly small in the z-axis, that is, the pertur- 
bation grows due mainly to the gas flow parallel 
to the slab. This tendency is more prominent 
when the slab is thinner. 



4.2. Non-Rotating Slab with Velocity 
Shear (f] = 0, [/ / 0) 

In the shearing slab both the gravitational 
and Kelvin-Helmholtz modes appear. The dis- 
persion relations described here are about the 
gravitational mode, which continuously approach 
those in the previous subsection 4.1, as the shear 
becomes weaker. (We will discuss briefly the 
Kelvin-Helmholtz modes in subsection 6.2.) 

Figure 5 shows the squared growth rate —uP' 
as functions of k^^ for three cases that Mn = 
uq{L)/c = 0,0.5, and 1.0, when L = 1.0 and 
ky = 0. Both — Wmax^ and fc^^max decrease as 
Mil increases. The dispersion relation remains 
parabolic-like on the oj"^ — k diagram. Figure 6 
shows the detailed behaviour of — Wmax^ against 
M|| for the gaseous slabs of L = 1.0, 0.5, and 0.3, 
when Q = and ky = 0. For a given L, —uJraax' 
is smaller when the shear is stronger. Particu- 
larly, when M\\ ^ 1, we have — Wmax^ ot Mm" 
(Compare with the thin line.). 

According to our extensive calculation the nu- 
merical dispersion relation can be well approxi- 
mated by 



Fig. 4. — The cross section of the slab perturbed 
by the fastest growing perturbation, when U = 



n = k 



y 



0. 
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(4.4) 



Equation (4.4) implies that the velocity shear 
may be dealt with as an ensemble of vortices and 
the gaseous slab has the effective sound speed, 

cJl + Mil 2 [See equation (4.1) 



Note that the 
square of the effective sound speed is the arith- 
metic sum of the square of the sound speed and 
that of velocity shear. From equation (4.4) we 



obtain 



and 
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Equation (4.5) and (4.6) are also consistent with 

-f^max^ Of M||~ when 



our numerical results, e.g. 
Ml, > 1. 



Fig. 5. — Dispersion relations when U ^ and 
n = 0. The thin line denotes the relation for the 
non-shearing slab. 



Figure 7 is the same as figure 4 except that the 
shear is present, kx = 0.285, ky = 0.0, L = 1.0, 
and Mil =1.0 are apphed. The concentration 
of gas is similar to that in figure 4, but the per- 
turbed gas spirals-in toward the point of maxi- 
mum density. Note, however, that the spirals are 
elongated along the z-axis. 

When r2 = 0, the wave number, {kx, ky) ap- 
pears in the matrix. A, and in the boundary con- 
ditions only in the form oik = Jkx'^ + ky'^ and 
kxUQ{z). Thus, we obtain a relation. 



ijj{kx, ky-, no, il = 0) 



Fig. 6. — The dependence of the squared maxi- 
mum growth rate on the shear magnitudes. The 



thin line is oj 



oc Mil 



uj{kx' = y . 



i^(f 1 rhy 5 i^y 



0;kxUo/kx',^ = 0) . 
(4.7) 

Using this relation we can obtain the growth 
rate for a mode having non- vanishing ky. Hence, 
the discussions about the dispersion relation (fig- 
ure 5) and the correlation between the maxi- 
mum growth rate and the Mach number of the 
streaming motion parallel to the slab (figure 
6) remain unchanged after the transformation 
from kx and Mm to k and Mukx/k. The effec- 
tive sound speed can be expressed as Ceff^ = 
c^ [1 + kx'^M\^/{kx^ + ky^)]. When fc^ = and 
ky 7^ 0, the growth rate is independent of the 
shear strength. 



Fig. 7. — The cross section of the shearing slab 
perturbed by the fastest growing perturbation. 



4.3. Rotating Slab without Shear (il / 

0,[/ = 0) 

Figure 8 shows the squared growth rates —oj"^ 
as a function of kx for ^2 = 0, 0.5, and 1.0. The 
parameters L and ky are kept constant, L = 1.0 
and ky = 0. The growth rate decreases mono- 
tonically, as 0, increases. The rotation tends to 
suppress the gravitational instability. However, 
the decrease in the growth rate is smaller than 
that for a rotating disk of which rotation axis is 
the normal to the disk. The latter growth rate is 
approximated by 



UJ 



2t 2 
C Kx 



2TTGakx + 417^ 



(4.8) 
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(see, e.g. Goldreich and TVemaine 1979) and the 
growth rate decreases by AQ,^ at a given kx- Ac- 
cordingly, as 0, increases, kxcr decreases while 
kxma.^ remains unchanged. On the other hand, 
kxraax decreases in our model as 0, increases. The 
critical wavenumber kxcr is independent of fi in 
our model. As proved in appendix, the velocity 
perturbation vanishes for the marginally stable 
mode of kx = kx cr in our model. Then the Cori- 
olis force does not work in the limit of iuj -^ +0. 
As 0, increases, the growth rate decreases but 
never vanishes as far as the mode is unstable at 
= 0. 

The correlation between maximum growth rate 
— i^max^ and the angular velocity, Q, of the slab 
is given in figure 9 for L = 0.1,0.3,0.5, and 1.0, 
when ky = and U = 0. As $7 increases, — Wmax^ 
decreases but a little and the decrease is very 
small for L < 0.1. 

These new results are intrinsic to the rotating 
slab sandwiched and pressed together with two 
plane-parallel rigid surfaces; when the gaseous 
component is perturbed, it moves easily along 
the slab, but the Coriolis force due to this motion 
drives the gas to move normal to the slab. The 
narrowly-separated rigid surfaces, however, sup- 
press this motion. The suppression is extremely 
strong, when the slab is thin, or the wavelength 
is long. Thus, the dispersion relations lack dra- 
matic effect of the rotation in that case. 

Figure 10a shows the perturbed density dis- 
tribution and is the same as figure 4 except for 
r^ = 0.4. We have taken a case of moder- 
ate slab thickness: L = 1.0, and assumed that 
ky = 0. Then, the maximum growth rate is 
~i^max^ = 0.25 and the corresponding wave num- 
ber is kxinax = 0.43. The global features are sim- 
ilar to those in figure 4, but we find that the 
iso-density contours are tilted due to the Corio- 
lis force. The velocity normal to the slab is ac- 
tually seen, derived by the Coriolis force. The 
concentration of gas is due primarily to the x- 
component of the perturbed velocity. 

As the angular velocity increases the iso-density 



Fig. 8. — Dispersion relations when 0,^0 and 
U = 0. The thin line corresponds to the non- 
rotating slab. 



Fig. 9. — The dependence of the squared maxi- 
mum growth rate on the angular velocity. 
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contours become more elongated and tilted, and 
then form two separated concentrations of gas 
(figure 10b). 

So far we have dealt with the case of kx 7^ 
and ky = 0. When ky ^ and k^ = 0, or 
the wave vector is parallel to the rotation axis, 
we find the growth rates hardly change due to 
the rotation, since the perturbed flow is nearly 
parallel to the rotation axis, and then 2Q x Wi ~ 
0. 



Fig. 10. — The cross section of the rotating slab 
disturbed by the fastest growing perturbation. 
The velocities seen from the rotating frame are 
delineated. We assume that fi = 0.4 and kx = 
0.43 in figure (a), and 0, = 1.0 and kx = 0.30 in 
figure (b). 



4.4. Rotating Slab with Shear {n^O,U ^ 
0) 

This case would be the most realistic gaseous 
slab occurring at the interface of the obliquely- 
colliding clouds with supersonic velocity. As seen 
from equation (3.13), the equilibrium state is 
valid only in the parameter range of QUH/c^ < 
1. Therefore, we present the results within it. 

Figure 11 shows the dispersion relation be- 
tween kx and — cj^ for Mm = 1.0, Q = 0.5, ky = 
0.0, and L = 1.0. For comparison, the dispersion 
relations for (M|| , n) = (0.0, 0.0), (1.0, 0.0), and 
(0.0, 0.5) are shown with thin curves for compar- 
ison. The wavenumber and the thickness of the 
slab are kept constant ky = and L = 1.0 for 
all the thin and thick curves. The gravitational 
instability is suppressed additively by shear and 
rotation, and the critical wave number kxcr de- 
creases as ft increases, which is derived to be con- 
stant when U = 0. 

Figure 12 shows the squared maximum growth 
rates as functions of fl for k 



y 



and L = 1.0. 

Three curves represent — Wmax^ for M|[ = 0, 0.5, 
and 1.0. In this figure we confirm again that 
the rotation and shear suppress the gravitational 
instability. Note that the maximum growth rate 
is larger for a larger Mm in the region of fi ^ 1.1. 
This is because the column density a is larger for 
a larger Mm for given L [See equation (3.15).]. A 
larger column density brings a larger maximum 
growth rate. 

Figure 13 is similar to both of figures 7 and 
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Fig. 11. — The dispersion relations when U ^ 
and Q ^ 0. Adopted Mm and Q for each hne are 
written in the figure. 



10a, but for My = 1.0, n = 0.5„ L = 1.0, k^ = 
0.2, and ky = 0. The contour curve of constant 
density slants to the x-direction as in figure 10 a 
and the gas stream spirals into the point of the 
maximum density as in figure 7. 

The dispersion relation between kx and —lo'^ 
in figure 12 cannot be empirically represented in 
a simple form including the parameters, the slab 
thickness 2L, shear U, and angular velocity fi. 
However, when the slab is a thin sheet of gas and 
the Coriolis force more suppressed, the dispersion 
relation reduced to that of the non-rotating slab 
in equation (4.4). 

5. Jeans Mass vi^ith Rotation and Shear 

When the gravitational instability grows and 
non-linear contraction proceeds, a high-density 
gaseous clump forms in the slab. When the 
slab rotates and has velocity shear, the Jeans 
wave length depends on the direction of the wave 
vector: The one parallel to the x-axis is given 
by Aa;j = 2Tr/kxmax and equation (4.6), when 
L <C -ffR, and the other one, A^j, parallel to the 
y-axis remain unaffected by the shear and rota- 
tion, and obtained from equation (4.3). There- 
fore, the most probable Jeans mass Mj would be 



Mj 



Aa;J Ay J 



cr. 



(5.1) 



Fig. 12. — The dependence of the squared maxi- 
mum growth rate on the angular velocity. 



where A^j = 27r//cj,max- Note that A^max, Ay^ax, 
and a are functions of the slab thickness 2L, 
which increases with time. Using a similar equa- 
tion, Whitworth et al. (1994) calculated the mass 
fragmenting in the shocked layer produced by the 
head-on collision. 

Using tm in equation (2.2), we define " the 
growth rate of the slab mass" as c^m = 1/tm- 
Assuming da/dt ~ a/t, where t denotes the time 
interval measured from the onset of the collision, 
we have 

(5.2) 



Wn 



1 



which monotonically decreases with time. We 
note, while, that the maximum growth rate. 
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— iwmax) of the instability increases with time, 
since a and L increase as the cohision proceeds 
[See also e.g. figure 3.]. At the beginning of the 
collision, uj^ » — ^"^max, and the slab mass grows 
more rapidly than the perturbation does, and 
thus the slab is practically stable. When — iwmax 
exceeds uj^-, the perturbation grows appreciably. 
Equating uj^ with — iwmax, we meet the time to 
when the perturbation starts to grow. 



*0 = l/l-it^max)- 



(5.3) 



We have then the Jeans mass Mj (5.1) in which 
^x:i-, Ay J, and a are evaluated at t = to- 

Using quantities defined in section 2, we have 



a 



M± ' 



(5.4) 



where we have used the relationships: a ~ 
tda/dt, da/dt ~ 2pocMj_, pi/po ^ M_i^, and 
Pi — Pci when L <C H^i. Substituting a in equa- 
tion (5.4) into (4.5), we have 



-ZUJ-n 



2ttGpc 



-.t. 



(5.5) 



M^^l + Mf 
From equation (5.3) and (5.5), to is written as 

1/4 



-(^)"'<-M„V' 



(5.6) 



Note that this is comparable to the free fall 
timescale, when M±_ ~ 1. Using the relations: 
kxraaa = vrG(T/c^(l + M||^) [from equation (4.6)], 
and kymux = vrGcr/c^ [from equation (4.3)], we 
have, respectively, 



cM_L 



(1 + M,| 



and 



Gpcto 
= 27rHM^^/\l + M\\^f\ (5.7) 

_ cM± _ 
Gpcto 

27^i7M_L^/^(l + M||2)~^/^ (5.8) 



where we have used equation (3.14). These are 
roughly comparable to or greater than H, when 
M± ~ 1. Finally, we obtain the Jeans mass from 
equation (5.1), 



Mj = 2tt^PcH^M±^/'^{1 + Ml 



2n3/4 



where we have used 



o-(t = to) = 2pcHM±-^/'^{l + Ml 



2^1/4 



(5.9) 



(5.10) 



Note that when M± ^ 1, we have L{t = to) <^ 
H ~ Hyi from equation (5.10) and a ~ 2pcL. 
When Mil tends to zero, equations (5.6) through 
(5.9) are reduced to equation (1) in Whitworth 
et al.(1994), except for numerical factors. 

The Jeans mass increases as U and fl increase, 
and therefore the more massive gaseous clump 
would be generated in the gaseous slab at the 
interface of the obliquely colliding gas clouds. 
When the collision is strong and the resultant 
gaseous slab is thin, the perturbed motion nor- 
mal to the shock front due to the Coriolis force is 
suppressed and thus the Jeans mass tends to that 
of the case of C/ 7^ and i7 = in subsection 4.2. 
We have derived it as analytical form in equa- 
tion (5.9). In this case, as is shown in equations 
(5.7) and (5.8), the clump formed by the gravi- 
tational instability is elongated along the x-axis 
and rotates around the y-axis. 

6. Conclusion and Discussion 

We have worked on the gravitational insta- 
bility of a gaseous slab with shear and rotation 
whose axis lies within and parallel to it. It is 
shown that more massive gaseous clumps are 
generated as the shear becomes stronger. The 
shear may be understood as ensemble of rotat- 
ing eddies, contributing to increasing the effec- 
tive sound speed. When the gaseous slab is 
thin and pressed together with two plane-parallel 
solid surfaces, the rotation affects little on the 
instability, because the perturbing force or the 
Coriolis force due to the rotation is directed only 
normal to them. We believe that the applied 
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plane-parallel solid surfaces are realistic enough 
to mimic the plane shock front which is stable 
against its deformation. We briefly touch upon 
below, what difference comes out from our results 
when we take the free boundary condition where 
the pressure always balance across the flexible 
and free surfaces. 



6.1. On the Boundary Condition Mim- 
icking the Shock Front 

In figures 2 and 3 we reproduce Elmegreen and 
Elemgreen's dispersion relation (1978) in thin 
lines, where [/ = fi = and the free bound- 
ary condition are adopted. The dispersion re- 
lations in figure 2 show that the difference be- 
tween the boundary conditions is more consid- 
erable as the slab thickness L decreases: com- 
pare (a,c) for L = 1.0 and (b,d) for L = 0.5. 
We find in figure 3 that the maximum growth 
rate — Wmax^ and their corresponding wavenum- 
ber are independent of the boundary conditions 
as far as L '^ 1. They, however, depend strongly 
on the boundary conditions when L <^ 1. When 
the boundary is "free " , — Wmax^ does not change 
in the wider range of L, while fe^max increases 
as L decreases (Elmegreen and Elmegreen 1978; 
Lubow and Pringle 1993). As many authors sug- 
gested, these properties seems due to the insta- 
bility whose density remain nearly constant while 
their configuration tend to round the structure 
(e.g. Vishniac 1983). Delineating the same fig- 
ure as figure 4 but now for the free boundary, we 
confirm this instability in figure 14. The veloc- 
ity field shows that the perturbed mass tend to 
round. 

The fixed and free boundary conditions would 
be two extreme cases representing a surface, since 
the former mimics completely rigid one, and the 
latter flexible one. The real structure of the 
shock front must be in between these two ex- 
tremes (See Welter 1982). We expect that the 
real boundary is closer to the rigid one when two 
clouds collide with high velocities and the shock 
wave is strong. 



Fig. 13. — The cross section of the shearing and 
rotating slab disturbed by the fastest growing 
perturbation. 



Fig. 14. — The same figure as figure 4 but for the 
free boundary. The wavenumber of this pertur- 
bation is now kx = 0.696 (and ky = d). 
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6.2. Kelvin-Helmholtz Instability in the 
Slab with Shear and No Rotation 

We have discussed exclusively the gravita- 
tional instability of the gaseous slab with shear 
and rotation. However, the basic equations in 
section 3 can automatically deal with the Kelvin- 
Helmholtz (K-H) instability, and it is apprecia- 
ble when the wave length of the perturbation is 
shorter, kx > 1, and the velocity gradient of 
the shear, duo/dz, is larger. In order to study 
the Kelvin-Helmholtz instability we made a slab 
model with strong velocity shear in which the 
unperturbed shear flow is assumed to be 



(6.1) 



uq{z) = [/tanh ( - 



M|| = 2.0, and n = 
is not measured in 

,2 



with / <C L. 

Figure 15 shows the relationship between kx 
and Re(-u;2), when ky = 0, L = 1.0, / = 1/30, 
0. We note that the ordinate 
-w^ as figure 2, 5, 8, and 11, 
but in Re (— w^), since —co"^ becomes the complex 
number. The solid lines denote the modes having 
real to"^ while the dashed lines denote the mode 
having complex w^. Thus each solid line is two- 
folded while each dashed line is four folded. In or- 
der to understand the network we labelled a solid 
line with a single character and a dashed line 
with two characters so that each symbol denotes 
a single continuous line. In the limit oi U ^> 
the zigzag line Ga — G — Gb — G — ■ approaches 
the dispersion relation of the gravitational in- 
stability. The sequences, Ga — a — ab — a — ■, 
■ — b — Gb — b — ab— ■, and ■ — c — ac — c — bc — c—- 
are the essentially the sound waves. These mode 
interact each other and form mixed modes, e.g., 
Ga, ab, be, etc. This mode coupling is common in 
the dispersion relation of the K-H instability in 
the slab with boundary surfaces (Glatzel 1987a, 
b). 

Figure 16 shows the density perturbations for 
the mode b of [kx, Re{-uj'^)] = (3.2, 0.62), the 
mode c of [k^, Re{-uj'^)] = (5.0, 0.47), and the 
mode 6c of [/c^, Ref-w^jj ^ (4^^ _i_q)_ We find 



Fig. 15. — Dispersion relations between kx and 
Re(— cj)^. In the segments delineated by thick 
solid lines, Imto = 0, while by thick dashed lines, 
Imoj 7^ 0. We calculate dispersion relations in a 
right region of the thin line. 



Fig. 16. — The cross section of the shearing 
slab disturbed by the growing sound waves. Iso- 
density contours are delineated. ky,L,l, Mu , and 
Q are the same as in figure 15. 
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that these have periodic structure along the z- 
direction with "wavelength" 2L/n : n = 2 for 
the mode h and n = 3 for the mode c. The 
mixed mode he have characters of both the par- 
ent modes, and is similar to the mode h in the 
region z > and while it is similar to the mode 
c in the region z < 0. 

Since the mixed mode has a complex fre- 
quency, they grow or decay while oscillating in 
time, i.e., propagating with the phase speed 
Rew/Zc^; in the x-direction. The density pertur- 
bation of the mixed mode has a smaller relative 
amplitude than the velocity perturbation of the 
mixed mode, |pi/po| < |i^i/c|. Thus, the growth 
of the mixed mode leads mainly to the diffusion 
of the velocity shear and does not likely to the 
formation of gravitationally bound systems. 

6.3. Non-Linear Evolution of the Unsta- 
ble Mass 

We have shown that the gravitationally unsta- 
ble clump collapses while rotating, when the slab 
has the shear or rotation. As many authors have 
pointed out, such a clump fragments into two 
parts during the collapse, and forms a binary sys- 
tem (or multiple system). Fujimoto and Kumai 
(1994) applies the results described in section 4 
for the formation of binary star-clusters observed 
in the Magellanic Clouds. 



Appendix. Proof of iti = at Critical 
Wave Number when iVf|| = 

The present appendix is to show that mi = 
at the critical wave number fccr, i-e., when the 
perturbation is marginally stable [uj = 0). We 
write the perturbed equations (3.16) through 
(3.20) in the non-dimensional form as in the text, 
and substitute equations (3.21) and (3.22) into 
them where we take a; = and wq = 0. We have 
then, writing k^r as h for simplicity. 



2^Uz\ = ikx —ipi , 

V PoJ 


(Al) 


= iky f-Vl - -) , 

V Po/ 


{A.2) 


d f pi\ 
2nuxi = -J- -V'l — , 
dz V PoJ 


{A.?,) 



ikxPoUxi + ikypouyi + -ripoUzi) = 0, {A A) 



and 



-^ - {kx + ky ) 



■l/jl = 2/91 



(^.5) 



[1) When ky ^ 0, we have from (A.2), 
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eral Scientific Research (04452013). 



Pi 



-poV'i- 



(^.6) 



Substituting equation (A. 6) into (A.l) and 
(A. 3), we have u^i = Uzi = since $7 / is pre- 
sumed. Using Uxi = Uzi = 0, we have Uyi = 0, 
because of /c„ 7^ in equation (A. 4), and hence. 



Ml =0. 



(^.7) 



(2) When ky = and kx 7^ 0, we combine 
equation (A.l) and (A. 3) to remove the term 
(-V'l - pi/po), 
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Uxl 



dz \ kx 



{A.i 



where we assume ^2 7^ as before. Substitut- 
ing equation (A. 8) into (A. 4), we have Uzi = 0. 
Then, equation (A.l) gives the relation (A. 6) 
again. Substituting equation (A. 6) into equa- 
tions (A.l) and (A. 3), we have 



Uxl = Uzi 



0. 



(^.9) 



Since we can choose Uyi = without loss of 
generality, we have again equation (A. 7). 
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